Product vacua with boundary states 
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We introduce a family of quantum spin chains with nearest-neighbor interactions that can serve to 
clarify and refine the classification of gapped quantum phases of such systems. The gapped ground 
states of these models can be described as a product vacuum with a finite number of particles 
bound to the edges. The numbers of particles, ul and ur, that can bind to the left and right edges 
of the finite chains serve as indices of the particular phase a model belongs to. All these ground 
states, which we call Product Vacua with Boundary States (PVBS) can be described as Matrix 
Product States (MPS). We present a curve of gapped Hamiltonians connecting the AKLT model 
to its representative PVBS model, which has indices ul = ur — 1. We also present examples with 
riL — ur — J, for any integer J > 1, that are related to a recently introduced class of SO(2J + 1)- 
invariant quantum spin chains. 



I. INTRODUCTION 

Gapped quantum phases and the transitions between 
them arc a subject of great current interest both for 
their fundamental importance, as illustrated by the ex- 
perimental observation of an _Eg-symmetry at the critical 
point of the quantum Ising model^, and because of their 
potential application in quantum information and com- 
putation as proposed in^. In the latter case the focus 
is on the topological phases^, which cannot be classified 
using a local order parameter, but depend on global fea- 
tures of the system such as the underlying topology of the 
model, as for example in the quantum Hall effect 4 . These 
topological orders can be characterized by an equivalent 
boundary theory^. 

Several recent works have aimed at classifying the 
ground state phases of quantum spin chains&i. While 
true topological order may not occur in one dimen- 
sion, direct analogues of some essential phenomena are 
present, including the role of edge states, which is the 
main topic of this article. 

The consensus is that the ground states of two mod- 
els Hq and Hi with short-range interactions are in the 
same phase if there exists a smooth path of Hamiltoni- 
ans H(s), < s < 1, H = H(Q),H 1 = H(l), such that 
the gap above the ground state does not close along this 
patb&£. Using the notion of automorphic equivalence, 
in£ we showed how this definition can be made precise so 
that it can also be applied for infinite systems. In general, 
proving that the gap does not close is a hard problem 
but a criterion exists for frustration free models whose 
ground states arc Matrix Product States (MPS)±£. Most 
progress has been made by focusing on such frustration 
free models, with or without prescribed symmetrie o 7 ' 11 . 

In this article, we introduce a new family of quantum 
spin chains with nearest-neighbor interactions with MPS 
ground states of a special type which we call Product 
Vacua with Boundary States (PVBS). The bulk phase is 
a product state (or several product states) , but for finite 
chains there are 2™ ground states that are obtained from 



the product vacuum by adding up to n distinguishable 
particles. Of these n particles til bind to the left edge of 
the interval and ur to the right edge. The equivalence 
classes of gapped ground states are in one-to-one cor- 
respondence with the values of the non-negative integers 
til and ur. The Hamiltonian preserves the particle num- 
ber for each type of particle separately. As an example, 
we identify the classes to which the SO(2 J + l)-invariant 
models ofi2 belong, namely Hl = ur = J. The AKLT 
model corresponds to J = 1. The complete class of PVBS 
models is more general and will be described elsewhere. 



II. PRODUCT VACUA WITH BOUNDARY 
STATES 

We consider a quantum spin chain with spin dimension 
d > n + 1. There are n + 1 states at each site that we 
interpret as n distinguishable particles and a vacuum. In 
addition there may be excited states of positive energy 
that will play no role in our discussion. For simplicity, 
we stick to the case d = n + 1. We start from a MPS 
representation of the ground states of the model. The 
Hamiltonian is constructed as the parent Hamiltonian for 
the set of MPS ground states and by general arguments 
we can then conclude that it has a gap and no other 
ground states^. 

Let label the empty state and let 1, . . . , n denote the 
n particle types. The ground states are generated by n + 
1 square matrices vq, v±, . . . , v n , satisfying the following 
commutation relations: 



ViVj = e 
v? = 0, i ^ 



2 



i ^ 3 



where 



'j 



and / A, € 



(1) 
(2) 
for 



< i , j < n . By redefining the phases 9ij, we can assume 
Xi > without loss of generality. We will also assume 
that Ao = 1, which amounts to a choice of normaliza- 
tion for vq. Although the exact form of the matrices 
is not essential for our discussion, it is important that 
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such matrices exists. One representation of the commu- 
tation relations is most conveniently described in terms 
of a chain of n spin- 1/2 particles, yielding matrices of 
dimension 2™. Define 

a+ = (o o) a = (? o) 
Wi = [o aJ P « = ( l] 

for z, j = 0, . . . n. The matrices are given by 

n 

Vq = P&Wi , 
i=l 

i—1 n 

Vi = PijWj ® cr + <g> PjfciUfe, i = l,...n. 

j=l fe=i+l 

Eq. |T]) follows from <7 + u>i = Ait«i(T + and PijU + = 
e l6i] /' 1 a + Pij as well as [Pjk,Wi] = for all k. Fi- 
nally, is a direct consequence of (c + ) 2 = 0. 

The MPS generated by this set of matrices for a chain 
of L spins are given by 

n 

1>{B)= n b(Bv iL ---v il )\i 1 ,...,i L ) , (3) 

»i,...,*t=0 

where B is an arbitrary 2™ x 2™ matrix. Consider the 
case L = 2. From the commutation relations (H} and @ 
it follows that all MPS vectors of the form ([3]) will be 
orthogonal to the vectors <pij G C d <E)C d given by 

<f> i = \0,i)-e- i9i °X i \i,0) (4) 

<ki = \h3)-e-" it K 1 *j\jrf (5) 
&i = |M) (6) 

for 1 < i,i < n and i j. If n > 2, not all combinations 
of particles can be realized in a ground state of a chain 
of length 2. Nevertheless, it is sufficient to consider a 
Hamiltonian with a nearest-neighbor interaction defined 
by 

n n 

z— i i<*<.7< n 

where • denotes normalization. The Hamiltonian for a 
finite chain of spins is then given by 

6-1 

H[a,b] = 2J ^.as+l) (?) 
a:— a 

where and h x ^ x+ \ is a copy of /i acting on the pair of spins 
at the sites x and x + 1. As a sum of orthogonal pro- 
jections, -ff[a,&] is non-negative and it is straightforward 
to verify that the MPS defined in ([3]) are eigenvectors 
with zero energy, hence ground states of the model. It is 
not hard to show that these are all the ground states of 

H[a,b] ■ 



If b—a+1 > n, for each subset . . . , i m } of {1, ... , n} 
there is a ground state fp^'^' 1 "" = tp(B %1 '---'' lm ) having ex- 
actly one particle of each type i\, . . . ,i m . The matrices 
B %1 '"' ,tm generating ^'••••'"» through (|3|) (up to normal- 
ization) can be chosen as 

/," ' =p ®(h-i) a - 01 ,®(i a -ii-i) ^-g, 

■ ■ • <g> <T _ <g> 

where p = cr+o' - . An interesting example is the ground 
state containing only particle i, 

b 

4, 6] = £(e ie <°A^|0,..., 0) (8) 

x—a 

where i is at site x in each term of the sum. We now 
add the assumption that Aj ^ 1, for 1 < i < n. The 
til particles having Ai < 1 arc bound to the left edge, 
whereas the tir = n — Ul particles with Ai > 1 are local- 
ized near the right edge, as can be seen in ([5]). The state 
with m = is the product state 

fi = |0, ...,0) . 

All other ground states differ from fi only near the edges. 
Specifically, 

6— >+oo 

for any local observable A. If only one of the edges is 
taken to infinity, the limiting ground states for the half- 
infinite chain depend on the particles at the other edge. 
Concretely, on the chain that extends to infinity on the 
right but with a left boundary, there remain 2 nL ground 
states corresponding to the possible combinations of the 
riL particles that bind to the left edge. Similarly, there 
are 2™ R ground states on the left infinite chain with a 
right boundary. 

Using the method oS^, we can prove that the energy 
of the first excited state is bounded below by a positive 
constant, independently of the length of the chain. As 
at most one particle of each type can bind to the edge, 
any second particle of that type must be in a scattering 
state. The dispersion relation of these scattering states 
can be explicitly calculated by considering the restric- 
tions of the Hamiltonian to any of the n invariant spaces 
containing exactly one particle. Properly centered and 
rescaled, these operators reduce to a free hopping Hamil- 
tonian with A-dcpcndcnt boundary conditions. A plane 
wave Ansatz yields the dispersion relation 

2A 

€i(fc) = 1 - ! 2 cos(fc + 9 i0 ) . 
1 + \ 

In particular the gap closes whenever Ai — > 1 for some 
i. Moreover, we conjecture that the exact gap in the 
thermodynamic limit is given by 

7 = nm t-r+AT :1 -'-"/■ 



3 



III. AUTOMORPHIC EQUIVALENCE AND 
GAPPED GROUND STATE PHASES 

Despite their simplicity, the PVBS models introduced 
in the previous section display the general characteristics 
of gapped one-dimensional systems with a unique bulk 
ground state and we can use them to illustrate the role 
of edge states in the classification of gapped phases. 

InJJi it is concluded that all gapped, translation in- 
variant, one-dimensional quantum spin systems without 
symmetry breaking belong to the same phase, and that 
they are equivalent to a product state. We believe that 
taking into account edge states, and in particular the 
behavior of the system on semi-infinite chains, is neces- 
sary as it allows for a finer classification closer in spirit 
to what one would find for models on two- and higher 
dimensional manifolds with nontrivial topology. 



A. The PVBS classes 

It is easily seen that two PVBS models of the type 
we introduced here belong to the same equivalence class 
if and only if they have the same values for the non- 
negative integers and tir. The reasoning is as fol- 
lows. Since equivalent phases are related by an automor- 
phism, a unique bulk ground state can only be mapped 
to another unique bulk state. Similarly, the ground state 
space dimensions of the half-infinite chains, 2" L and 2 nR , 
are also preserved by an automorphism. Hence, if two 
PVBS models belong to the same phase, they must have 
equal nj, and tir. Conversely, if two PVBS models have 
the same values of Ul and tir but each with their own 
sets of parameters {\i(a) | I < i < til + tir} and 
{Oij(a) | I < i,j < til + tir}, for a = 0,1, an inter- 
polating path along which the gap does not close can be 
constructed as follows. First, one may apply a strictly 
local unitary to perform a change of basis in spin space 
such that both are PVBS expressed in the same spin basis 
and such that A^ (a) < 1 for 1 < i < nj, and Xi (a) > 1 for 
hl + 1 < i < til+ur, for both a — and a — 1. Let u be 
the unitary for this change of basis. Then, take a smooth 
curve of unitaries u(s), < s < 1, with u(0) = 1 and 
u(l) = u, and let J7[ a ,6](s) be the (6 — a + l)-fold tensor 
product of u(s). Now we conjugate the initial H[ a b ] with 
these unitaries to define a smooth path of Hamiltonians 
with a constant gap. Simultaneously, we can deform the 
parameters of the two models by linear interpolation: 

Ai(s) = (l-s)A i (0) + sA i (l), 
%( S ) = (l-s)%(0)+ S %(l). 

This yields a smooth family of vectors 4>ij( s ) as m (@][6|) 
and thereby a smooth family of nearest-neighbor inter- 
actions h(s) and of Hamiltonians. The gap remains open 
because A;(s) ^ 1 for all 1 < i < n and s e [0, 1]. By 
the general result of£ this implies the quasi-local auto- 
morphic equivalence of the two models. Note that it is 
essential that in each pair between which we interpolate, 



the Ai's are either both < 1 or both > 1, which is why 
we had to assume that and tir are the same for both 
models. If one uses the same type of interpolation to con- 
nect models with different values of and hr, the gap 
necessarily closes along the path and there is a quantum 
phase transition. This is not to say that one could not 
construct paths along which the gap closes also for the 
case of constant and tir and of course this would then 
not imply a transition between different gapped phases. 

The uncountable family of PVBS models which depend 
on the real parameters {A^} and {%} is completely clas- 
sified by the pair of integers (tiljUr). Given their sim- 
plicity, it is natural to choose them as representatives of 
the much larger phase they belong to. 



B. The AKLT model 

As a first example, we show that the AKLT models 
belongs to the same equivalence class as the PVBS mod- 
els with til = tir = 1. The AKLT model is an antifcr- 
romagnetic spin-1 chain with a unique, gapped ground 
state in the thermodynamic limit, and four zero-energy 
states on a finite chain, which are usually described in 
terms of a spin- 1/2 particle attached to the two ends of 
the chain. Wc found a smooth path of gapped mod- 
els connecting the AKLT model with a PBVS model 
with one particle for each boundary. Let us denote the 
two particle states by — and +. For s £ [0, sq] where 
sin(so) = i/2/3, the following four vectors span the 
ground state space of two neighboring spins of the the 
interpolating models as a function of s: 



^ 00 (s) =n(s)sm(s) [A(f 



V>- + (s) 



-cos 2 ( S )(l + A( s ) 4 )|0,0) 
-A( S )|0,-) + |- 0) 
-A(s)|+,0) + |0,+) 
|-,+)-A(.s) 2 |+,-), 



where A(s) is a smooth function such that A(so) = 1, < 
A(s) < 1, for all s < s , and n(s) = (1- A(s) 2 cos 2 (s)) 1/2 . 
The corresponding nearest-neighbor interaction, 



h(s) = i ]T 

(u)G{(00),(0-),(0+),(-+)} 

is the projection onto the orthogonal complement of 
this 4-dimcnsional space. We then define H[ ab ](s) — 

X^=i hx.x+i(s)- It readily follows that -H] a ,&]( s o) is the 
AKLT Hamiltonian and that H [afi] (0) is the PVBS model 
with til = tir = 1, the coefficients A_ = A(0) and 
A+ = A(0) _1 , and all the phases 9^ — tt. The path 
of interactions is smooth as the four ground state vectors 
are smooth, remain orthogonal to each other and of fi- 
nite norm for all s, and the spectral gap does not closed. 
Hence, the AKLT model belongs to the same gapped 
quantum phase as the PVBS model with = ur = 1. 
In particular, the sets of ground states of these models 
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are automorphically equivalent for the finite, half-infinite 
and infinite chains, where they are isomorphic to a pair of 
qubits, a single qubit, and a unique pure product state, 
respectively. 

The ground states of -ff[ a ,b] ( s ) f° r s > have a minimal 
matrix product representation with 2x2 matrices Vi(s) 
defined as follows: 




By contrast, the MPS representation of the ground states 
of the PVBS model, at s = 0, uses 4x4 matrices. In 
fact, there is no faithful 2-dimcnsional representation of 
the commutation relations (JTJ [2]) with n = 2. Up to 
unitary change of basis, the only nilpotent matrix is er + , 
so we can set v+ = c+(7 + and V- = c~U*a + U , with 
c± =/= and U =/= 1. But their commutation relation im- 
plies Tv(v-V + ) = so that U = 1 if c± ^ 0. The parent 
Hamiltonian constructed from the limiting matrices Vi (0) 
corresponds to a PVBS model with til = 1 arid tir = 0. 
With a minor modification one can interchange the roles 
of til and Tin, but in either case the limiting matrices 
lead to a parent Hamiltonian that is not in the same 
phase as the AKLT model; It turns out to be equivalent 
only in the bulk and on the half-infinite chain tending to 
oo on the right but not on the chain that is infinite on the 
left. We believe that the finer classification we propose in 
this paper yields interesting additional information com- 
pared to previous approaches such as**. We conclude 
that it is essential to construct a smooth, gapped family 
of Hamiltonians, and that it is not sufficient to find a 
smooth family of MPS matrices only. A similar conclu- 
sion was reached iniS, where 'uncle Hamiltonians' were 
constructed that share the ground state spaces of gapped 
'parent Hamiltonians' but are gapless in the thermody- 
namic limit. 

C. The SO(2J + 1) models 

The S'J7(2)-symmetric spin-1 AKLT chain has been 
generalized in a number of different directions: higher 
S'C/(2)-spinsi^, higher lattice dimensions^, SU(N)- 
invariant models with N > 321, and most recently to 
a class of spin-J models with an SO(2J + l)-invariant 
nearest-neighbor interaction^. The simplest way to in- 
troduce the latter generalization is to note that the ker- 
nel of the interaction of the AKLT model, i.e., the ground 
state space of two neighboring spins is spanned by the an- 
tisymmetric vectors, namely the spin I triplet, and spin 
singlet state, which is symmetric. As shown inl£, The 
antisymmetric subspace of two spins enhanced with one, 
arbitrary, symmetric state, is the ground state space of 
a frustration-free spin chain with a translation invariant 
nearest-neighbor interaction. In particular, the spin-J 
chain, with integral J > 1, with the nearest- neighbor 



interaction given by the projection onto the orthogonal 
complement of the span of the antisymmetric vectors and 
the (symmetric) spin singlet state is a frustration-free 
model with a unique gapped MPS ground state of the 
infinite chain. It is easy to see that this interaction com- 
mutes with SO(2J + 1) acting by its fundamental repre- 
sentation on each spin. 

In the case of the AKLT model (J = I), the antisym- 
metry of the ground states is reflected in the fact that the 
matrices (j^lllOp at s = so can be related to the generators 
of the Clifford algebra C3 : 

Z 1 = y/3/2(v+ - v-), Z 2 = - V372i(w+ + «-) (11) 
Z a = -\f3v , 

In fact, the transformation {v-,vq,v+} — > {Zq,Zi,Z 2 } 
corresponds, up to a globl rescaling by v3, to a change 
of basis in the local Hilbcrt space of the spins. Simi- 
larly for J > 1, the ground states of the SO (2 J + 1) 
model are, up to an overall normalization factor, gener- 
ated by the matrices {Z a = Z* a : < a < 2 J} of the 
higher dimensional Clifford algebra C 2 j+i satisfying the 
anticommutation relations 

ZaZp + ZpZ a = 28 a @l. (12) 

The symmetry of the ground state is manifest in the 
invariance of (|12p under the transformation Z'^ = 
^2 a Op a Z a for any orthogonal matrix O. The normal- 
ization factor 7 can be found by setting ^2 a j 2 Z^ = 1, 
yielding 7 = (2 J + l)" 1 / 2 . 

A representation of the Clifford algebra C2J+1 can be 
obtained from a representation of the algebra of canon- 
ical anticommutation relations (CAR) with J creation 
operators 




where q = 1 - p = a a + . Note that (— l) q = a 3 . Then, 
for 1 < j < J, 

Z-y-i = a j + a * 3 7 Z 2 j = i(ctj — a*) , 

and 

J J 
Z = Y[(2a*a 3 -l) = (g)(-iy. 

3=1 3=1 

The canonical anticommutation relations of the 
(a*,aj)j =1 imply the Clifford relations for the {Z a )^L . 
This representation is of dimension 2 J . 
Now, the matrices 

Vzj-i = ajcij V 2 j = Pja* V = 7Z0 

arc related to the normalized Clifford generators by a 
change of basis if and only if \a 3 ; | 2 = \/3j | 2 = 2/ (2 J + 1). 
In that case, the set {V } U {V 2 j-i,V 2 j : 1 < j < J] 
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generate the same matrix product states as the (Z a )^L , 
namely the ground states of the SO(2J + 1) invariant 

model. Note that JV I a j| 2 + l7| 2 = lj so that it is natural 
to set 



7 =: cos(sq), and 



aij = -^=sin(s ). 



for all j = 1, ... , J. 

It remains to introduce a deformation of the CAR alge- 



bra to relate it to the PVBS algebra with nj, = 
We first note that for any complex number A, 



A • a~ X q = X q c 



a+X q = A • AV 



where X q = p + Xq. For parameters Ai, . . . , Xj. 
duce the following twisted creation operators 



;(A) = (-A 1 )«g).. 

their adjoints, and 
do (A) = 



(-A 



'A,V 



-n R = J. 

, (13) 
we intro- 



(-Ai) 5 



Clearly, a* (A) 2 = = ay (A) 2 . A direct consequence 
of (Till)) arc the twisted commutation relations: 

a*(A)a,(A) + A~« , ; A :■</ ; •: A ; = a (A) 2 
a* (AH (A) + Xj X k a k (X)a* (A) = 



0' * *) 



a*(A)a (A) + A J a (A)a*(A) = 
a*(X)al(X)+X j Xfa* k (X)a*(X)=0 



and their adjoint relations. The higher dimensional ana- 
log of the algebraic path (JSJ [TU)) is easily obtained with 
the following definitions for s G [0, So]: 

Vjy-x(s) = aj(s)aj(s), V 2j (s) = /3j(s)a*(s) 

and Vq(s) = j(s)a (s), where a!j(s) = a*(A(s)) for a 
smooth path of parameter vector A(s) such that Xj(s) < 1 
for s < s and Aj(so) = 1. Using the commutation rela- 
tions, the normalization ^ Q y a Ta = 1 reads 



i = E(i&i 2 -A?k-i 2 ) %-«i+[Ek-i 



which implies 



J 

E 

3=1 



a/ + | 7 | 2 = l 



and 



=1-A 2 (1- 



Note that the SO(2J + 1) symmetry is broken as soon 
as |j9j| 7^ | ctj ] . Concretely we choose 7 (s) = cos(s) and 
Oj(s) = sm(s)/y/j for s € [0, So], thereby producing a 
path that mimics (J9J 1 10[) . The commutation relations 
of the V a 's interpolate between those of the CAR (and 
therefore of the SO(2J + 1) models) and of the PVBS 
with parameters 



In particular, 

V 2j (s)V 2j -i(s) + XjisfVv-iWVijis) 



sin( S ) [l-A 2 (l-sin 2 ( S )/J)] T; 



yfj COS 2 (s) 



V s 
v o ■ 



{Aj(0), A^O) -1 : 1 < j < J}; 



(14) 



This algebraic path of matrices generates a smooth 
path of matrix product states, from which a smooth 
family of Hamiltonians with nearest-neighbor interaction 
can be constructed, for any J. They interpolate be- 
tween the SO(2J + 1) model and the PVBS model with 
iil = riR = J. The properties needed to conclude that 
the gap remains open along the path can be derived from 
a suitable generalization of the arguments for the AKLT 
path, J = 1, and will be included in future work. 



IV. DISCUSSION 

The fact that one of each of the n types of particles 
can appear in a ground state of the PVBS system, i.e., 
without raising the energy, should not be interpreted as 
implying that the particles are massless. Quite to the 
contrary, there is a mass gap for each of them. It turns 
out, however, that they can bind to the left or right edge 
of the chain and the binding energy exactly equals the 
mass gap so that such states with a particle bound to 
the edge are degenerate with the vacuum ground state. 

As we explained in the previous section, the AKLT 
model is equivalent to a PVBS model with two types of 
particles, one that binds to the left edge and one that 
binds to the right edge. For a finite chain, this yields 
a four-dimensional ground state space (no particles; one 
particle on the left boundary; one particle on the right 
boundary; two particles, one on the left and one on the 
right). The bulk ground state is unique. It has a fi- 
nite correlation length in the AKLT model but is a sim- 
ple product state in the equivalent PVBS model. The 
unitary transformation relating the two models is quasi- 
local^, meaning that, under conjugation with the unitary, 
local observables map to observables that are approxi- 
mately local, i.e., that up to an arbitrarily small correc- 
tion depend only on a finite number of spins. Since our 
goal is to understand the structure of bulk phases and 
edge states (and, in higher dimensions also the effects 
of topology), the quasi- locality of this unitary transfor- 
mation is essential. In finite volume, any unitary will 
preserve the dimension of the ground state space. The 
importance of the locality property comes into sharp fo- 
cus when one takes the thermodynamic limit and consid- 
ers the resulting ground states on the infinite and half- 
infinite chains. For example, the non-local unitary trans- 
formation introduced by Kennedy and Tasaki 2 ^ to reveal 
the hidden string order— transforms the 4 AKLT ground 
states into 4 translation invariant product states, which 
leads to 4 distinct bulk ground states in the thermody- 
namic limit. This example shows that non-local unitary 
transformations do not preserve the structure of the bulk 
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ground state(s) and clearly would not be useful to classify 
gapped ground state phases. 

The bulk ground state of a SO(2J + 1) model for any 
integer J is equivalent to a unique product state. This 
however is not sufficient for them to all belong to the 
same gapped phase as the number of edge states depends 
on J. The definition of a phase through a quasi-local 
automorphism yields a finer classification that takes these 
boundary states into account. 

Several generalizations of the PVBS are possible. 
First, it is straightforward to describe models with more 
than one vacuum state, e.g., models with a broken dis- 
crete symmetry. This could include breaking of the trans- 
lation symmetry of the lattice to a subgroup, leading to 
periodic ground states. Second, the construction can eas- 
ily be generalized to allow more than one particle to oc- 
cupy the same site. The auxiliary chain of n spin-1/2 
particles we have used in the appendix just has to be 
replaced by an arbitrary finite spin chain. Third, if all 



or a subset of the A-parameters are equal, one can de- 
scribe Hamiltonians with a local continuous symmetry 
group. Adding symmetry as a constraint considerably 
enriches the classification problem of the gapped ground 
state phases and the phase transitions between themii^i. 

It appears to us that the class of models with PVBS 
ground states is able to capture the range of behaviors 
seen in gapped ground states of one-dimensional spin sys- 
tems with short range interactions. The most interest- 
ing phenomena from the physical point of view occur in 
two dimensions. Some aspects of the classification prob- 
lem have already been considered in two and higher di- 
mensions. The notion of automorphic equivalence stud- 
ied ins, e.g., applies to arbitrary dimensions. Interest- 
ing boundary states in two dimensions have been studied 
ir*^. Much further work is needed to better understand 
the landscape of gapped ground state phases in two and 
three dimensions. 
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